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SUMMARY 

A  method  is  developed  for  the  minimization  of  gear  tooth  surface  deviations 
between  theoretical  and  real  surfaces  to  improve  the  precision  of  surface  manufac¬ 
ture.  Coordinate  measurement  machinery  is  used  to  determine  a  grid  of  surface 
coordinates.  Theoretical  calculationa  are  made  for  the  grid  points.  A  least-square 
method  is  used  to  minimize  the  deviations  between  real  and  theoretical  surfaces  by 
altering  the  manufacturing  machine-tool  settings.  An  example  is  given  for  a  hypoid 
gear . 


INTRODUCTION 

The  Gleason  Works  were  pioneers  in  the  application  of  coordinate  measurements  to 
improve  the  precision  manufacturing  of  hypoid  and  spiral  bevel  gears  (ref.  1).  In 
aerospace  applications,  duplication  of  flight-qualified  master  gears  is  very  impor¬ 
tant,  and  coordinate  measurement  has  now  become  part  of  the  normal  production  pro¬ 
cess.  Methods  to  enhance  and  extend  the  use  of  this  machinery  can  be  very  valuable 
to  aerospace  gear  manufacturers. 

The  approach  developed  in  this  paper  enables  one  do  determine  deviations  of  a 
real  surface  from  the  known  theoretical  surface.  This  is  accomplished  by  using 
coordinate  measurements  and  minimizing  the  deviations  to  correct  the  previously 
applied  machine-tool  settings.  The  surface  deviations  are  represented  in  the  direc¬ 
tion  of  the  normal  to  the  theoretical  surface.  The  coordinate  measurements  are 
performed  by  a  machine  with  4  or  5  degrees  of  freedom.  In  the  case  of  4  degrees  of 
freedom,  the  probe  performs  three  translational  motions  (fig.  1);  the  fourth  motion, 
rotation,  is  performed  by  turning  the  table  with  the  workpiece.  The  axis  of 
rotational  motion  coincides  with  the  axis  of  the  workpiece.  In  the  case  of  a 
5-degree-of-f reedom  machine,  the  fifth  degree  of  freedom  is  used  to  provide  the  probe 
deflections  in  the  direction  of  the  normal  to  the  theoretical  surface.  The  probe  is 
provided  with  a  changeable  spherical  surface  whose  diameter  can  be  chosen  from  a  wide 
range . 

The  motions  of  the  probe  and  the  workpiece  for  coordinate  measurements  are  com¬ 
puter  controlled  and  for  this  purpose  a  grid,  the  set  of  points  on  the  surface  to  be 
measured,  must  be  chosen.  There  is  a  reference  point,  one  point  on  the  grid,  that  is 
necessary  for  the  initial  installments  of  the  probe.  There  are  two  orientations  of 
the  probe  installment  that  are  applied  for  measurements  of  a  gear  (fig.  1(a))  and  a 
pinion  (fig.  1(b)),  depending  on  the  anglt  £  iue  p...Lv,ii  cone. 


Probe 


Gear-\ 


(a) 


Figure  1 — Surface  measurement  orientation  for  (a)  gear 
and  (b)  pinion. 

The  mathematical  aspects  of  coordinate  measurements  will  now  be  described 
(ref.  2):  First,  it  is  necessary  to  derive  the  equations  for  the  theoretical  sur¬ 
face.  In  many  cases,  this  surface  can  be  derived  as  the  envelope  to  the  family  of 
generating  surfaces,  namely  the  tool  surfaces.  Next,  the  results  of  coordinate 
measurements  must  be  transformed  into  deviations  of  the  real  surface  represented  in 
the  direction  of  the  surface  normal.  Here,  the  surface  variations  are  represented  in 
terms  of  the  corrections  to  the  machine-tool  settings.  The  surface  deviations 
obtained  from  coordinate  measurements  and  the  surface  variations  determined  by  the 
corrections  of  machine-tool  settings  can  be  represented  by  an  overdetermined  system 
of  linear  equations.  The  number  k  of  these  equations  is  equal  to  the  number  of 
grid  points,  and  the  number  of  unknowns  m  is  equal  to  the  number  of  corrections  of 
machine-tool  settings  (m  «  k) .  The  optimal  solution  to  such  a  system  of  linear 
equations  enables  one  to  determine  the  sought-for  corrections  of  machine-tool 
settings. 


EQUATIONS  OF  THEORETICAL  TOOTH  SURFACE 

Considering  that  the  theoretical  surface  can  be  determined  directly,  we  repre¬ 
sent  it  in  coordinate  system  ’n  two-parametric  form  as 

r^{u,9),  n^(u,0)  (1) 

where  and  are  the  position  vector  and  unit  normal  to  the  surface,  respec¬ 

tively,  and  (u,0)  are  the  Gaussian  coordinates  (surface  coordinates). 


2 


For  the  case  when  surface 
surface  E  ,  we  represent  in  S 
(ref.  3) 


is  the  envelope  to  the  family  of  generating 


surface 


and  the  unit  normal 


to 


as 


(2) 


=  [Lt^]nc(u„,e^),f{u^,0c/^)  =0  (3) 

where  (u  ,6  )  are  the  Gaussian  coordinates  of  the  generating  surface  E  ,  and  ^  is 
the  generalized  parameter  of  motion  in  the  process  fox"  generation.  The  equation  of 
meshing  is 


«  , a  Jk . 

'  ''C  '  ''  C  '  ^  ' 


ret) 

V 


=  0 


where  is  the  normal  to  E  ,  and  is  the  relative  motion  for  a  point  of 

contact  of  E  and  E  .  Matrices  (M  1  and  [L  1  describe  the  coordinate 
transformation  from  S  to  for  a  position  vector  and  surface  normal,  respec¬ 

tively.  Position  vectors  in  three-dimensional  space  are  represented  with  homogeneous 
coordinates . 


COORDINATE  SYSTEMS  USED  FOR  COORDINATE  MEASUREMENTS 

Coordinate  systems  and  S^  are  rigidly  connected  to  the  coordinate  measur¬ 

ing  machine  (CMM)  and  the  workpiece  being  measured,  respectively  (fig.  2).  The 


Figure  2  —Relationship  between  theoretical  and  measurement 
coordinate  systems,  (p,  radial  distance  to  point  from  axis  of 
rotation.) 
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Figi  're  3 — Measurement  grid  on  tooth  surface. 


backface  cf  the  gear  is  installed  flush  with  the  base  plane  of  CMM.  The  distance  1 
between  the  origins  O  and  O  is  assumed  to  be  known,  but  the  parameter  6  of 
orientation  must  be  determined  (see  the  following  section).  The  coordinate  transfer- 
mation  from  to  is  represented  by  the  matrix  equation 

(5) 

m  ^  mt  t 


MEASUREMENT  GRID  AND  ESTABLISHMENT  OF  THE  REFERENCE  POINT 

The  grid  is  a  set  of  points  on  chosen  as  points  of  contact  between  the 

probe  and  (fig-  3).  Fixing  the  value  of  for  the  point  of  the  grid  and  the 

value  of,  say,  y^  (or  obtain  the  following  equations; 

y,(u.,0J  =  h.,  z^(Uj,0.)  =1.  (i  =  l,...,k)  (^) 

where  k  is  the  number  of  grid  points. 

We  consider  h^  and  1^  as  given  and  solve  equations  (6)  for  (u^,0^).  Then  we 
can  determine  the  position  vectors  and  the  unit  normals  for  k  points  of  the  grid 
using  the  equations 

r!'’  =  [x^(u^,0j)y^(u^,0.)z^(u^0.) ]^,  (i=l,...,k) 


(i  =  l,...,k) 

The  position  vector  for  the  center  of  the  probe,  if  the  deviations  are  zero,  is 
represented  by 
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(i)  (i) 


(i) 


(i  =  1, . . . ,k) 


(9) 


where  p  is  the  radius  of  the  probe  sphere. 
The  reference  point 


(o) 


=  [x^(u  ,d  )yt(u  )z^(u  ) 


(10) 


is  usually  chosen  as  the  mean  point  of  the  grid. 

The  center  of  the  probe  that  corresponds  to  the  reference  point  on  is 

determined  from  equation  (9)  as 


(o) 


o'"'  rv 


(o)  nO)  , 


(o)  n(o) ,  T 


(11) 


where  (u*°’,6'°’)  are  known  values. 

The  coordinates  of  the  reference  center  of  the  probe  are  represented  in  coordi¬ 
nate  system  S  of  the  measuring  machine  by  the  matrix  equation 


(o)  (O) 


(12) 


Equation  (12)  yields 


(o)  (O)  c  (O)  n(o) 


<°l  ,  c  (o)  n(o)  . 

y  =  y  (u,u  ,u  ) 


(O)  (O)  c  (O)  Mo) 

z  =  z  ( o , u  ,u  ) 


(13) 


-  fo)  (o)  (o) 

Three  equations  (13)  contain  four  unknowns:  v,  '^111  '  solve  these  equa¬ 

tions,  we  may  consider  that  one  of  the  coordinates  of  the  reference  point  of  the 

(  O) 

probe  center,  say,  y^  ,  may  be  chosen  equal  to  zero.  This  is  accomplished  by 
requiring  the  reference  point  to  lie  in  the  x  -  z  plane.  The  orientation  of 
angle  0  is  now  established  to  satisfy  this  requirement,  and  all  measurements  are 
referenced  from  this  location.  Then  equation  system  (13)  allows  one  to  determine 

J-  (o)  .  ^  ,  (o>  (°)  „  (°) 

0,  X  and  z  (ref.  2).  Coordinates  x  ,y  =  0,  z  are  necessary  for  the 
m  m  m  -'m  in  ^ 

initial  installment  of  the  center  of  the  probe. 
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MEASUREMENT  OF  THE  DEVIATIONS  OF  THE  REAL  SURFACE 

The  deviations  of  the  real  surface  are  caused  by  manufacturing  errors,  heat 
treatment,  etc.  Vector  positions  of  the  center  of  the  probe  for  the  theoretical 
surface  and  the  real  surface  can  be  represented  aa  follows: 


R  =  r  (u,0)  +  Xn  (u,0) 

where  r  and  n  are  the  position  vector  and  the  unit  normal  to  the  theoretical 

mm* 

surface  and  are  represented  in  coordinate  system  S  of  the  measuring  machine;  A 

in 

determines  the  real  location  of  the  probe  center  and  is  considered  along  the  normal 

♦ 

to  the  theoretical  surface;  R  and  R  represent  in  S  the  position  vector  of  the 
probe  center  for  the  theoretical  and  real  surfaces,  respectively.  Equations  (14) 
and  (15)  yield 

R  -  R  =  (X  -  p)n  =  Ann  (1^) 

m  m  '  r  '  m  n\ 


An  =  (R  -  R  ) 


The  position  vector  R^  is  determined  by  coordinate  measurements  for  points  of  the 
grid.  Equation  (17)  determines  numerically  the  function 

An^=Anj^(u^,0^)  (i=l,...,)c)  (18) 

that  represents  the  deviations  of  the  real  surface  for  each  point  of  the  grid. 


MACHINE  TOOL  SETTINGS  TO  MINIMIZE  DEVIATIONS 


The  procedure  used  to  minimize  the  deviations  can  be  represented  in  two  stages: 
(1)  determination  of  variations  of  theoretical  surface  caused  by  changes  of  applied 
machine-tool  settings,  and  (2)  minimization  of  deviations  of  real  surface  by  appro¬ 
priate  correction  of  machine-tool  settings. 

We  consider  that  the  theoretical  surface  is  represented  in  S^  as 

r^  =  r^(u,0,dp  ( j  =  1, . . . ,m)  (19) 

where  parameters  d  are  the  machine-tool  settings.  The  surface  variation  is 
represented  by 


5r,  =  _(5u  *  ^  66  *  T  _(5d, 

"  du  W  “  ad^  ^ 
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We  multiply  both  sides  of  equation  (20)  by  the  surface  unit  normal  and  take  into 

account  that  5r  /50  •  n  =  9r  /9u  -0=0,  since  9r  /90  and  9r  /9u  lie  in  the 

t  1-  t  >■  t  t 

plane  that  is  tangent  to  the  surface.  Then  we  obtain 


n  =  T)  -  *  Ilf  5d  =  a5d 


We  can  now  consider  a  system  of  k  linear  equations  in  m  unknowns  (m  <<  k)  of  the 
following  structure: 


a,,6d,  +  a,.,(5d,  +  .  .  .  +  a,  6d  =  b, 

11  1  12  2  lin  in  1 


a^,(5d,  +  a,,5d,  +  .  .  .  +  a.  5d  =  b^ 

kl  1  k2  2  Km  :r  k 


(22) 


where  i  designates  the  number  of  grid  points;  a  .  (s  =  l,...,k;  j  =  l,...,m) 
represents  the  dot  product  of  partial  derivatives  9r^/9d,  and  unit  normal  n  .  The 
system  (22)  of  linear  equations  is  overdetermined  since  m  <<  k.  The  essence  of  the 
procedure  for  minimization  of  deviations  is  determining  unknowns  5d.  (j  =  l,...,m) 
that  will  minimize  the  difference  between  the  left  and  right  sides  of^equations  (22). 
The  solution  employed  the  least-square  method.  The  subroutine  DLSQRR  of  IMSL  MATH/ 
LIBRART  (ref.  4)  was  used  to  computerize  the  procedure. 


APPLICATION  OF  METHOD  TO  THE  INSPECTION  OF  FORMATE  HYPOID  GEAR 

Each  tooth  side  of  a  formate  face-hobbed  gear  is  generated  by  a  cone,  and  the 
gear  tooth  surface  is  the  surface  of  the  generating  cone.  Two  cones  that  are  shown 
in  figure  4(a)  represent  both  sides  of  the  gear  space.  The  following  equations 
represent  in  coordinate  system  S  gear  surfaces  for  both  sides  and  the  unit  normal 
to  such  surfaces  (fig.  4(b)); 

-Sgcos 

(r  -  SgSin  ag)sin  6^ 

(r  -  SgSin  aQ)cos  9^ 

1 
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sin  a.. 


-  cos  a„sin  8^ 

O  -j 

-  cos  a,,COS  6^ 


where,  r  is  the  position  vector  and  the  unit  normal;  r  is  the  cutter  tip 

radius;  Ct,  is  the  cutter  blade  angle  (3^  >  0  for  the  concave  side  and  a^,  <  0  for 
the  convex  side). 

Figure  5  shows  the  installment  of  the  generating  cone  on  the  cutting  machine. 
Coordinate  systems  and  are  rigidly  connected  to  the  cutting  machine  and  the 

gear  being  generated,  respectively.  Systems  S  ,  S  ,  and  are  rigidly  connected 


to  each  other  since 
gear  tooth  surface 
equat ions : 

the 

gear  is  formate  cut 
and  the  unit  normal 

To  represent  in 
to  E_,  v;c  use  the 

the  theoretical 
follov.'ing  matrix 

r,(Sg,0^,dj  = 

[M,jr^(s^,0^) 

(26) 

n^(s^,0,.,d,)  = 

[L,^]n^(s,.,0,,) 

(27) 

where 

[M.J  = 


cos  7^ 

0 

-sin  7^ 

1  m 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

-V, 

sin  7 

f  W 

0 

oos  7„ 

E 

X 

<3 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

The  surface  Gaussian  coordinates  are  s^.  and  8^,  and  d^(7^^,  V^,  H^,  and  ^te 

the  machine-tool  settings. 

The  numerical  example  presented  in  this  paper  is  based  on  the  experrmein.  that 
has  been  performed  at  the  Dana  Corporation  (Fort  Wayne,  IN,  U.S.A.).  The  initial 
deviations  An  for  each  side  of  the  real  tooth  surface  have  been  obtained  by 
measurements  on  a  coordinate  measuring  machine  (fig.  6).  The  grid'  for  the  meas¬ 
urement  is  formed  by  nine  sections  along  the  tooth  length,  each  section  having  five 
points.  The  number  of  grid  points  k  is  therefore  45,  and  the  reference  point  is  at 
the  middle  of  the  grid,  i.e.,  the  third  point  of  the  fifth  section.  In  the  measure- 

(  O  ) 

ment,  the  coordinate  y^  of  the  reference  point  is  chosen  to  be  zero  and  the 
alignment  angle  6  is  determined  from  solving  eguation  system  (13). 

The  minimization  of  deviations  was  performed  in  accordance  with  the  algorithm 
described  in  MACHINE  TOOL  SETINGS  TO  MINIMIZE  DEVIATIONS,  and  the  results  are 
illustrated  in  figure  7  and  table  T. 
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Figure  7  Minimized  deviations  after  corrections  made  to  machine-tool  settings 


90 


TABLE  I.  -  RESULTS  OF  MINIMIZATION 

[Pressure  angle,  a,  =  21.25°;  cutter  diameters  =  9  in.; 
point  width  of  cutters  =  0.08  in.] 


Machine-tool 
sett ings 

Machine-tool  setting  parameters 

1 

mm 

mm 

1 

ry 

rad 

Ax  , 

mm 

Initial 

Corrected 

103 .252550 

103 .25220 

27.466600 

27.21603 

1.059816 

1.06437 

0.009677 

-0.53343 

CONCLUSION 

A  general  approach  for  a  computerized  determination  of  deviations  of  a  real  sur¬ 
face  from  the  theoretical  one  based  on  coordinate  measurements  has  been  proposed.  An 
algorithm  for  computerized  minimization  of  deviations  by  corrections  of  initially 
applied  machine-tool  settings  has  been  developed.  The  approach  is  illustrated  with 
the  example  of  the  tooth  surface  of  a  hypoid  formate  gear. 
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